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DKfilGN OF MDR MODULATORS 



The present invention relates to r^harmacophores 
which are useful for identifying and designing molecules that 
possess mult i- drug resistance (MDR) modulatory activity. 

The clinical resistance of tumours to treatment with 
certain cytotoxic agents is an obstacle to the successful 
chemotherapeutic treatment of cancer patients. A tumour may 
acquire resistance to a cytotoxic agent used in a previous 
tv-eattnent. A tumour may also manifest intrinsic resistance, 
or cross-resistance, to a cytotoxic agent to which it has not 
previously been exposed, that agent being unrelated by 
structure or mechanism of action to any agent used in 
previous treatments of the tumour. 

Analogously, certain pathogens may acquire 
resistance to pharmaceutical agents used in previous 
treatments of the diseases or disorders to which those 
pathogens give rise. Pathogens may also manifest intrinsic 
"resistance, or cross resistance, to pharmaceutical agents to 
which they have not previously been exposed. Examples of 
this effect include multi-drug resistant forms of -anlnr.^ , 
tuberculosis, leishmaniasis and amoebic dysentery. 

The above phenomena are referred to collectively as 
multi-drug resistance (MDR) . A 170KDa plasma membrane 
glycoprotein (P-gp) is implicated in the mechanism which 

underlies MDR. 

Cancer cells which exhibit multi-drug resistance, 

referred to as MDR cells, display a reduction in 
intracellular drug accumulation compared with the 
corresponding drug-sensitive cells. Studies using in vitro 
derived MDR cell lines have shown that MDR is often 
associated with increased expression of ?-gp which has drug 
binding properties. P-gp is thought to function as an efflux 
pump for many hydrophobic compounds, and transfection studies 
using cloned P-gp have shown that its overexpression can 
confer the MDR phenotype on cells; see, for example, Ann. 
Rev. Biochem 58 137-171 (1989). 

A major function of P-gp in normal tissues is to 
export intracellular toxins from the cell. There is evidence 
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to suggest that overexpression of P-gp may play a clinical 
ro^e in multi-drug resistance. Increased levels of P-gp mRNA 

o-otein have been detected in many forms of human cancers 
- leukaemias, lyrr.phomas, sarcomas and carcinonas. Indeed, m 
sor'.e cases P-go levels have been found to be increased m 
tumour bioosios obtained after relapse from chemotherapy. 

inhibition of P-gp function in P-gp mediated MDR has 
been shown to lead to a net accumulation of anti-cancer agent 
the cells. For example. Verapamil a knovm calcium channel 
block'^- v;as shown to sensitise MDR cells to Vinca alkaloids 
in vi^ro and in vivo: Cancer Res., 41, 1967-1972 (1981) . The 
oroposed rr.echanisn of act.on of the MDR modulatory compound 
Involves ccu.F'^c .taon with the anti-cancer agent for bxnd.xng 

to the P-gp- ^ 

The ;m modulatory compounds can therefore be used 

in a T-.echod of potentiating the cytotoxicity of an agent 

which is cvtccoxvc to a tumour cell. Such a method 

comnrises.'fcr instance, administering an MDR modulatory _ 

comoound to .h. :u-our cell whilst tho tumour c.ll is exposec 

to 'Che cvtoto:<;c ..gent in question. The therapeutic effect 

of a ch.^^rch=..-acc.utic. or antineoplastic, agent may thus be 

e-ranc^d. T.. ■ r.-.l t . -drug resistance of a tumour ceil to a 

cytotoxic .-.vor.: uurxng cherr.otherapy may be reduced or 

eliminated. u ^ 

r/cd-yiatorv comoounds can also be used in a method of 
treating a disease in which the pathogen concerned exhibits 
multi-drug resistance, for instance multi-drug resistant 
forms of malaria IPlasjr^dimD ialciEor^ tuberculosis, 
leishmaniasis and a^.ccbic dysentery. Such a method 
comprises, for instance, administering an MDR modulatory 
compound with (separately, simultaneously or sequentially) 
the drug to which the pathogen concerned exhibits multi-drug 
resistance. The therapeutic effect of the drug may thus be 

enhanced. . 

A human or animal patient harbouring a tumour may be 

treated for resistance to a chemotherapeutic agent by a 

method comprising the administration thereto of an MDR 

modulatory compound. The MDR modulatory compound is 

administered in an amount effective to potentiate the 



cytotoxicity of the said chemotherapeutic agent. Examples of 
chemotherapeutic or antineoplastic agents which are preferred 
in the context of the present invention include Vinca 
alkaloids such as vincristine and vinblastine; anthracycl ine 
antibiotics such as daunorubicin and doxorubicin; 
tnitoxantrone; actinomycin D; taxanes e.g. taxol; 
epipodophyllotoxins e.g. etoposide and plicamycin. 

In addition, a human or animal patient suffering 
from a disease in which the responsible pathogen exhibits 
mult i -drug resistance may be treated for resistance to a 
therapeutic agent by a method comprising the administration 
thereto of an MDR modulatory compound. 

Examples of such disease include multi-drug 
resistant forms of malaria ( Plasmodium falciparum) , 
tuberculosis, leishmaniasis and amoebic dysentery. 

MDR modulators also have utility in the delivery of 
drugs across the blood-brain barrier, and in the treatment of 
AIDS and AIDS-related complex. A human or animal patient in 
need of such treatment may be treated by a method comprising 
the administration thereto of an MDR modulatory compound. 

Disadvantages of drugs which have been used to 
modulate MDR, termed resistance modifying agents or Ri-lAs, 
include that they frequently possess a poor pharmacokinetic 
profile and/or are toxic at the concentrations required for 
MDR modulation. 

Molecules that possess MDR modulatory activity span 
a diverse range of structural types. 

The elucidation of which functional group or 
collection of functional groups bestows a certain biological 
activity on a small organic molecule has long been c^ 
interest to the pharmaceutical industry. Historically, 
however, the lack of satisfactory answering has relegated the 
development of active compounds to good fortune or exhaustive 
synthesis and biological testing of a large number of 
compounds. With the advent of computer aided drug design 
this situation is beginning to change. As an increasing 
amount of information is being gained about drug receptors 
and the conformations that are energetically accessible to a 
drug molecule at normal (physiological) temperature, a more 



complete understanding is being developed of the drug - 
receptor interaction. 

It is well understood that, on binding, molecules 
recognise each other by the molecular potenti.>l fields that 
they generate as they approach one-another. The field is not 
only dependenr. on the molecular structure but also on the 
molecule's conformation on approach to the receptor molecule 
and the induced field changes that the approach causes. 
Thus, it has been suggested that field correspondence rather 
than structural correspondence should be a guide to 
similarity betv;een active molecules. 

Attempts to define molecular fields in terms of 
electrostatic parameters have been made, for example, using 
CoMFA and QSAR techniques. These are restricted in the 
number of conformations they can handle, are bound to fixed 
grids and ignore interstices or use energy calculations 
parameterisGd for isolated molecules. They are 
compatutionally demanding and can not generally account for a 
diverse range of structures having the same bioactivity. 
Thus, using such methods, a reliable pharmacophore has not 
yet been developed, particularly not for MDR activity. 

The invention provides a pharmacophore coiviprisin^a 
3 -dimensional array of field points defining a shape and 
volume •.vheroin the set of 3 -dimensional field points is the 
aggregate average of the field points derived from a 
plurality of molecules which possess MDR modulatory activity. 

The iiiventioxT also provides a method of screening a 
molecule for MDR modulatory activity, which method comprises; 

(i) obtaining field points describing the positive, 
negative and van der Waals surface around the molecule; 

(ii) measuring the energy of overlay when the molecular 
field points obtained in (i) are fitted into a pharmacophore 
of the invention; and 

(iii) selecting the molecule as a candidate MDR modulator 
if the energy of overlay measured in (ii) is equal to or 
lower than the energy of overlay between field points for 
molecules known to possess MDR modulatory activity and the 

said pharmacophore. 

The invention also provides a method of designing a 



molecule having MDR modulatory activity, which method 
comprising the step of determining the structure of a 
compound which has field points that have, when fitted into a 
pharmacophore of the invention, an energy of overlay equal to 
or lower than a predetermined value. 

The invention also provides a compound which 
possesses MDR modulatory activity identified using a 
pharmacophore of the invention or by a method of the 
invention, or designed using a method of the invention. 

A pharmacophore can be thought of in general terms 
as any 3 dimensional array of descriptors which go towards 
explaining the therapeutic activity of a drug or receptor. 

A pharmacophore as used herein, is defined as a 
three dimensional array of field points defining a shape and 
volume. The field points are indicative of areas of positive 
charge, areas of negative charge and areas of van der Waals 
interactions that are present in the field points derived 
from active compounds. 

By positive field point is meant a point or area in 
space around a molecule which has electrophilic properties. 
A positive field point is most likely to interact with a 
negative field point on another molecule, for instance on 
binding to a receptor protein. 

By negative field point is meant a point or area in 
space around a molecule which has nucleophilic properties. A 
negative field point is most likely to interact with a 
positive field point on another molecule, for instance on 
binding to a receptor protein. 

By van der Waals surface field point is meant a - 
lipophilic region that has a strong possibility of close 
approach to the receptor macromolecule in this case P-Cp. By 
close approach it is meant a distance such that van der Waals 
forces are able to come into effect. 

Fig 1. Depicts a conformation of XR9051 with XED charges 
added . 

Fig 2. Depicts a conformation of XR9051 with related field 
points. Field colours being: 
Red « positive, 
green « negative and 



yellow --^ van der Waals regions. 
Fig 3. Shows a conformational mapping diagram describing 

the top 3 links between molecules. 
Fig 4. Depicts the optimum cyclic path and related 

conf ormers . 

Fig 5. Depicts a representation of a typical pharmacophore 

of the invention. 
Fig 6. Shows a plot of log [accumulation IC^o] against 

coulonibic overlay energy for a number of knovm 

modulatory compounds. 

Pharraacophores have previously been defined as three 
dimensional arrangements of chemical groups related to a 
aiven biological activity which enables meaningful comparison 
of molecules exhibiting the same biological function (Naruto 
5 et al. Eur. J. Med Chem'. 20:529.532 (1985)). Pharmacophores 
thus defined can be derived by simple overlap of active 
structures or common functional groups in the molecules. 
However, a problem is encountered when molecules which 
c:uiibic a great diversity of structures and structural mota f s 
0 are oroved to be active for the same biological function. 

Simply overlapping the active structures or common functional 
group's in the m.oiecuies does no- produce a reliable 
phari>iacophore. Simple pharmacophores such as tliese 
frequently do not take into account the conformation of. or 
15 energetic accessibility to, the 'active' molecule. Without a 
way to orient the molecules e.g., based upon fit v;ith another 
macromolecule - a receptor, enzyme, or in this case 
P-gp it is difficult, and in some cases impossible, to 
construct a reliable pharrr.acophore . This problem results, at 
20 least in part, from the fact that sometimes even closely 
related active molecules fit into macromolec ^les m very 

different ways. 

In the present invention a pharmacophore is a 

three dimensional array of field points. Therefore a 

25 pharmacophore may be thought of as representing an aggregate 

array of positions in space of the field points of a series 

of molecules having the same or similar biological activity. 

Pharmacophores can be specific for different compounds, their 

related molecules and a particular biological activity. In 



the present invention the pharmacophore is specific for MDR 
activity. It should be emphasized that these created 
pharmacophores do not exist as such in nature and are the 
product of aligning the field points of a plurality of 
molecules . 

A pharmacophore, once created, stands alone and is 
subsequently independent from the molecules, in the present 
case MDR m.odulatory molecules, that were involved in its 
formation. Thus, after formation of the pharmacophore, one 
no longer needs to use structural comparison with the 
molecules themselves for the design or identification of 
biologically active molecules. A pharmacophore of the 
invention can be used to identify and design new molecules 
that will possess the same or similar structural and charge 
features that are represented by the pharmacophore. This is 
a completely different concept from the one of using the 
simple structural similarities between active molecules as 
the\asis for the design of compounds. The pharmacophore can 
be Vised itself for any number of applications, including but 
not limited to the following: as a screening tool for drug 
develoomont and to determine if a particular compound will 
possess bioactivity of a certain type, in the present case 

lADR modulatory activity. 

Each pharmacophore has characteristic constraints 
placed upon its shape, topology, volume, and electrostatic 
profile. A Dharmacophore is accurately described by its 
three dim.ensional array of field points which is represented 
by a coordinate system that is configured in a computer 
memorv. Each specific atom within a molecule that fits in a 
> pharmacophore has a specific location relative to the active 
site of the macromolecule with wh^ ch it is to interact. The 
individual atoms also have electrical fields assigned to 
them. These fields are represented numerically and through 
many other ways including the use of colours and shading to 
5 indicate the field strength. As the degree of steric and 
electrostatic fit between the pharmacophore and the field 
points derived from a candidate molecule, whether said 
candidate molecule has previously been synthesised or not, 
increases, a more negative coulombic overlay energy results. 
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This indicates an increase in the efficacy of the 
pharmacophore for that candidate molecule. This could 
manifest itself as increased bioactivity. 

The term "coulombic overlay energy" as used herein 
is representative of the total energy of fitting between two 
field patterns, each field pattern being a field point 
descriptor file. Thus the greater the negative value of 
coulombic overlay energy, the better the fit of the two field 
patterns. A successful pharmacophore would be expected to 
produce a good correlation between the coulombic overlay 
energy when a molecule is fitted into the pharmacophore and 
the MDR modulatory activity of the molecule. 

The pharmacophores of the present invention can be 
generated using any commercial or research software capable 
of the processes defined below. Pharmacophores of the 
invention are preferably generated using the "cosmic" 
molecular modelling suite comprising "xedminlin" , 
"xedf Idprep" , "xedminOl", "xedqmfcom", "proc_bmf", "xedsort" 
and "xeddraw" running on a computer suitable for running 
molecular modelling software being equipped with appropriate 
accessories. Preferably Silicon Graphics Indy. R4400, 
150MHZ. 64MB RAM, 2GB HD equipped with XZ Graphics running 
IRIX 5.3 operating system. 

The "cosmic" modelling suite is copyrighted to 
Dr J G Vinter, Department of Chemistry. University of 
Cambridge, Lensfield Road, Cambridge, England, CB2 lEW. 

Pharmacophores of the present invention are derived 
from the structures of molecules which are known to 
possess MDR modulatory activity preferably molecules which 
are known to possess potent MDR modulatory activity, more 
preferably molecules known to inhibit multi-drug resistance 
in vivo and/or in vitro and are more preferably still 
molecules including at least one of XR9051, M3-073. GF120918, 
CP100356. S9788. VX-710 and BB509-035. 

The 3-dimensional array of field points of which the 
pharmacophores of the invention are comprised are derived 
from the structures of the pluarity of molecules which 
possess MDR modulatory activity preferably by: 
(i) subjecting each molecule to conformational hunting 



to generate a plurality of conformers for each molecule; 

(ii) applying extended electron density (XED) charges to 
each molecular conformer generated in (i) ; 

(iii) defining field points describing the positive, 

5 negative and van der Waals surface around tach conformer; 

(iv) making an intermolecular energy comparison between 
the field points defined in (iii) ; 

(V) obtaining the three lowest energy comparisons for 

each intermolecular comparison and mapping them on a 
10 conformational mapping diagram; 

(vi) defining an optimum cyclic path which relates the 

molecules; and 

(vii) defining the field points for the conformer of each 
molecule which lies on the optimum cyclic path defined in 

15 (vi) . 

The structures of the plurality of molecules known 
to possess MDR modulatory activity are provided via 
construction with a chemical structure drawing program 
preferably compatible directly or indirectly with the 
20 torsional minimization and field preparation program 

"xedminlin". Preferably the structures are provided via the 
program "xeddraw" as cosmic.dat files. Each molecule is 
independently subjected to conformational hunting preferably 
using the program "xedminlin" . 
25 If the program "xedminlin" is used then: 

the bond torsion increment angle is preferably 5- 

360«; more preferably 10«>; 

carbon to carbon double bonds are preferably set as 

non rotatable bonds; 
30 amide carbon to nitrogen bonds are preferably set as 

trans and as non rotatable; 

if a flexible ring structure is present then the 
number of randomisations of that ring is preferably 5; 

the upper limit of energy above the global minimum 

35 is preferably 20 Kcal; 

the number of randomisations of the whole molecule 

is preferably 200; and 

the dielectric constant is preferably 1 to 4, more 

preferably 1. 
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The results of this process are a set of energy 
minimized conformations for each molecule. If using 
-.xedminlin" then these are presented in a cosmic. bst file. 

Extended electron density (XED) charges are then 
5 added to each of the conformers of each of the molecules 
preferably using the program "xedf Idprep" which is called 
automatically by "xedminlin". Further information about the 
calculation of XED charges from first principles can be found 
in J.G. Vinter. J . Comput . -Aided Mol. Design 8 (1994) 653-668 
10 and J.G. Vinter et al, J. Comput . -Aided Mol.Design 9 (1995) 
297-307 and references cited therein. 

Field points are then defined for each of the 
confonnations of each of the molecules preferably using the 

program "xedf Idprep" . 
15 If "xedfldprep" is used to generate field pomto for 

each conformation then the probe used in field generation is 
preferably set to "peptide /peptide interface". 

The information defining each field point comprises 
the type of field point; positive field point, negative field 
point and van der Waals surface field point; the 3D 
co-ordinates defining the position of the field point and a 
numerical value quantifying the size/extent of the field 
point . 

A number of lowest energy conformations for eacti 
25 molecule are selected for intermolecular comparisons. The 
number of conformations is preferably selected in the 
following way: 

(i) Selecting the 40 lowest energy conformations tor 

each molecule, then 
30 (ii) selecting only the compounds from (i) which have an 

energy of 8 Kcal above the global minimum (lowest energy 
conformer) for the molecule then 

(iii) reducing the upper energy limit to 3 Kcal above the 

global minimum, and then, - , . „ 

(iv) further raising or reducing the upper energy limit 
so that a convenient number of conformers are selected. The 
number should be chosen so that the intermolecular comparison 
that follows can be carried out in a practical amount of 
time . 
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35 
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Most preferably the number of conformations selected 
for each molecule is from 5 to 25 e.g. from 7 to 21. 

If using "xedfldprep" then the conformation with XED 
changes and field points added is generated as a cosmxc.fst 

Comoarisons betv;een pairs of molecules are then made 
solely on the basis of the field point descriptors (type, 
location, extent) for each conformation of each of ^r.e 
molecules. The comparisons are preferably made using the 
0 method of least squares to establish the best fit between tne 
^ield ooints of each of the two conformers. The comparisons 
Ir. preferably made using the program "xedqmf com', which is 
called automatically on selecting the BCOM menu item from .h. 
-xedfldprep^' program. When using ••■xedcn.f com" comparisons are 
5 ouantified by coulomioic overlay energy. '^^^ ^^^^^^^^^^//^ 
of one set of field points over the other. The better 
the f^t the lo-.ver the coulom:oic overlay energy. For eacn 
inter^olecular comparison the 50 having the lowest coulombic 

overlay energy are selected. .^ori a 

■U using "xcdqn,£conV then the data is generated as a 



20 

cooiaic.bmf file 



ihe results of the comparisons are then processed. 
orefe-ably ^^;-3ing the program "proc_bmf". If u.3ing ^'^roc.hn-.f ' 
;>.e ^o.ult of each intermolecular comparison between eacn 
25 conformer is displayed as coulombic overlay -ergy. 

I- using "oroc_bmf" to process the comparison data 
then a text file displaying the comparisons is generated; a 

• fiiP is created automatically 

cosmic.com. Also an ^.quivalent file is crea 

which is compatible with the next stage of pharmacophore 

dPGTun; a cosmic, sot file. 

In.arcon£or.aticnal comparisons between conformers 
of pa.rs of molecules are ranked in order o£ coulombic 

, .. The results o£ the comparisons 

overlay energy, lowest first. Ine result , „„„ 

are .napped onto a conformational mapping diagram * 
35 sLple conformational napping diag:.am can be obtained if only 
^he top three ranking comparisons are mapped. Each conformer 
selected for use in the intermolecular comparisons is 
assigned a nu^er. lumbers representing the conformers whicn 
are in the top three ranked comparisons for each pair of 
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compounds are placed in a pattern or shape e.g. around a 
circle. Pairs of conformers which are in the top three 
ranked comparisons are joined, each line representing a 
coulon^bic overlay energy for the comparisons of that pair of 
conformers. (see figure 3). The conformational mapping is 
carried out using a suitable mapping program, preferably the 
prograra "xedsort". When using the program "xedsort" it is 
preferable to select the top 3 ranked comparisons between 

each molecule. 

The optimum cyclic path between molecules which 
relates specific conformations is then defined using a 
suitable computer program, preferably "xedsort" . 

A cyclic path is a path which can be followed on the 
conformational mapping diagram between all of the molecules 

5 represented >-:oarc from one to another wherein one 

conformation of .^ach molecule lies on the path. The optimum 
cvclic pach iu th3 cyclic path which has the lowest sum of 
coulombic cv3rUv/ energies. The optimum cyclic path is found 
automatically if -.xodsort" is used. In this case all cyclic 

0 paths U-..-C- .ionr.-I .-r.d thoir total coulombic overlay energies 
calculated. Tr.j lowest is then selected (see figure 4). 

Th- ccntormacions that fail on this optimum cyclic 
path are up:.r.-c ch-- -active" conformations of each of the 
molecules. 

5 Tr. ' fio.c i;oints for each conf or;r,ation which lies on 

the optimur- cycUr path are thus known. If using the cosmic 
modellinr: sui^c then each conformation is located in a 

cosmic. fsc flic. 

A phai-scophore is then generated by combining the 
0 field ooin-.s of the" ^'active" conformations using an aggregate 
average method. The parameters, such as the area over which 
a field point is averaged, can be adjusted to define 
differing pharmacophores. The pharmacophore is generated 
using a suitable computer program, preferably the "create 
>5 template" feature of the "xedsort" program. When using the 
"xedsort" program for generating the pharmacophore tne 
preferred parameters are: 

Electronic Hunting Distance (EHD) = 2.0 to 3.0 
surface Hunting Distance (SHD) = 0.5 to 3.0 
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SHD 



= 2.0 to 
3 . 0 and 



5.0 



Lower Limit Filter (LLF) 
most preferably EHD =3.0, 

LLF =2.0 

5 The pharmacophore is generated ac a cosmic. cpf file. 

The most preferred pharmacophore is a pharmacophore having 
the following field points: 
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15 



20 



25 



30 



35 



Atom 
Type 

-7 

-7 

-7 

-7 

-7 

-7 

-7 

-7 

-7 

-6 

-6 

-6 
-6 
-6 
-6 
-6 
-6 
-6 
-6 
-6 
-6 
-6 
-6 
-6 
"-6 

--6 
-6 
-6 
-6 
-6 



•1. 

4 . 
-4 

3. 

7 
-6 



X 

381791 
,879077 
.418141 
. 562676 
.206078 
,492228 
-6.899206 
-4 .444944 
-7.304530 
-6.831137 
5.201669 
4 .797647 
4 .724869 
-1 .074833 
5.790484 
-5.261753 
-1 .473675 
-5.881620 
0.611350 
-5.783878 
1.734998 
7.598676 
9.587200 
8.219660 
4.675870 
-1.133006 
-7.544412 
-9.125527 
-7.119754 
2.332763 



Y 

-0.051772 
-2.537416 
0.089392 
4 .774324 
1.839680 
-10.028677 
-6.415226 
-7.302401 
8.956919 
-4 .965901 
0.095870 
5.131311 
0.054068 
-1.256638 
3.896064 
1.922220 
5.651587 
-0.674191 
-0.417927 
3.835954 
-0.534123 
-2.714864 
-2.109952 
0.966784 
2.821225 
-9.683890 
-1.565759 
-5.337989 
-8.783302 
-5.381765 



Z 

1.715311 
-5.803063 
6.975651 
6.500367 
- 1 . 274007 
1.675134 
2.683271 
-3.535873 
0. 177705 
5. 561902 
-2.055756 
-2 . 954147 
5.399575 
-0 .461180 
1.509006 
4 .363908 
0.238804 
-5.583631 
5.020914 
-4 .612012 

10.273318 
-1.632125 
-4.796533 
-8.134950 
-7.802368 
0.147326 
0.572997 
-1.040235 
-1.247938 
3.370511 



Charge 

-14.7214 
-6.0644 
-4.3289 
-8.6588 
-8. 5662 
-2.1906 
-7 .7529 
-4.6301 
-2.2923 
2.3235 
8.2589 
6.4282 
5.1285 
8 .2910 
3.2358 
6.6584 
6.5755 

3.9685 
7.3589 

6.1148 
2.7267 
4.56^- 
4.9354 
5.0419 
4.9785 
5.2926 
3.8282 
2.6320 
3.5238 
12.2559 
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-6 -6.407646 

-6 -1.476288 

-6 -9.308277 

-6 -2.701459 

-6 -0.104383 

-5 0.780038 

-5 -6.232273 

-5 -3.467500 

-5 5.578010 

-5 -0.161253 

-5 -2.238659 

-5 -0.744474 

-5 -1.279436 

-5 -7.480459 

-5 -1.934318 

-5 -3.708762 

-5 3.115965 

-5 7.338803 

-5 2.050096 

-5 --1.193728 
-4 .348258 
-0 .521176 



■2 . 
-4 , 
4 

-6 



,912221 
.227778 
.502877 
.564950 
5.439674 
0 .691121 
2 .601393 
-2 .270888 
-4 .081225 



-3.392133 
1.803737 
0.694358 
8 .367678 
7 . 873184 
1 , 882956 
1.696347 
3 . 913391 
0.238672 



0.037560 -6.019390 
-0.779217 8.042641 



-6 .033126 
-5.982479 
-2 .520825 
6.335472 
5 .991115 
-2.167775 
2 .696229 
4 .896177 



-0. 571065 
3 . 824223 
0 . 141700 
6.394043 
-0.484472 
-0 . 311654 
-2 . 834882 
5. 076291 



-5 
-5 



3 .6672 

7 .4663 
2 .8107 

8 .4013 
2 .2347 
-11 .5525 
-10.7761 

-11.2323 
-3 .5105 
-8.1238 
-8 . 1531 

-8.8733 
-10 .0235 

-8.1691 
-5.5024 
-6 . 8736 
-3 . 9490 
-5 .4513 
-7 . 9992 
-7 . 2592 
-8 . 1012 



-2.180720 -3.354927 
1.087337 -5.092674 
2.639100 10.414149 -3.5695 
wherein acor.^. cype describes the nature of the field poi:.t, 
-7 represents surface field points, -6 represents positive 
field poincs and -5 represents negative field points; X, Y 
and Z represent co-ordinate values; and charge represents the 

effective size of the field. 

The phanracophores of this invention can be sho'.-n to 
possess a good relation with MDR modulatory activity. This 
car. uc demonstrated by comparing the pharmacophore with the 
field points of each of the plurality of compounds from whxch 
the Dharmacophore is derived (already available as a 
cosmic. f St file) using the overlaying comparison method usea 
previously. The comparison of the field points of each 
cofapound with the pharmacophore gives a coulombic overlay 
energy value. This value can then be correlated with the 
bioactivity of the compound. The most preferred way of 
measuring how well the pharmacophore correlates to activity 
is to plot log [accumulation IC„1 against coulombic overlay 
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calculating the sample correlation coefficient (r) • 
Preferred pharmacophores of this invention, when exposed to 
this procedure, have an (r) value of not less than 0.5 
preferably not less than 0.65. 
5 Pharmacophores of this invention can be used as a 

tool for screening a molecule for MDR modulatory activity. 
One such method of screening a molecule for MDR modulatory 
activity comprises; 

(i) obtaining field points describing the positive, 
10 negative and van de Waals surface around the molecule; 

(ii) measuring the energy of overlay when the molecular 
field points obtained in (i) are fitted into a pharmacophore 
of the invention; and 

(iii) selecting the molecule as a candidate MDR modulator 
15 if the energy of overlay mentioned in (ii) is equal to or 

lower than the energy of overlay between field points for 
molecules known to possess MDR modulatory activity and the 
said pharracophore . 

Pharmacophores of this invention can also be used as 

20 a tool for designing a molecule having MDR modulatory 

activity. One method of designing a molecule having MDR 
modulatory activity comprises the step of determining the 
structure of a ;jmpound which has field points that have, 
when fitted into a pharmacophore of the invention, an energy 

25 of overlay equal to or lower than a pre-determined value. 

The pre-determined value can be arbitrarily set or can be the 
energy of overlay measured when the field points for a 
plurality of molecules known to possess MDR modulatory 
activity are fitted into the said pharmacophore. 

3 0 Compounds which possess MDR modulatory activity can 

be identified using a pharmacophore of the invention. MDR 
modulatory compounds can be identified by the method 
mentioned herein or by any other method using pharmacophores 
of this invention. Compounds which possess MDR modulatory 

35 activity can be designed either by using a method as 

mentioned herein or by any other method using pharmacophores 
of the invention. 



.'J 
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Examples 

The invention is illustrated by way of the follov;ing 
examples with reference to Figures 1 to 6 - 

The pharmacophores were produced using "cosmic" 
5 molecular modelling suite comprising "xedminlin", 

"xedf Idprep" , "xedminOl" , "xedqmfcom" , "proc-bmf " , "xedsorc" 
and "xeddraw", running on a Silicon Graphics Indy, R44 00, 
150MHZ, 64MB RAM, 2GB HD equipped with XZ Graphics and 
running IRIX 5.3 operating system. 
10 Example 1 

The following potent MDR modulators which exhibit a 
wide range of structural types were used to produce a 
oharmacoohore of the invention. 
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The structure of each molecule was independently 
drawn using "xeddraw" . Each molecule was then subjected to 
conformational hunting using "xedminlin" . The starting 
structures were imported into "xedminlin" as cosmic.dat 
files. The parameters were set as follows :- 

Bond increment angle 

Upper energy limit above global min 

Number of randomisations 

Dielectric constant 

Number of randomisations of flexible ring 
Hold carbon carbon double bond as entered Yes 
Hold amide carbon nitrogen bond as entered Yes 
Probe for field generation: peptide/peptide 

interface . 

The program "xedminlin" automatically conformation 
hunted the molecules, added XED charges and calculated the 
field ooints for 'the 40 lowest energy conf orm.ations . A 
number" of lowest energy conformations were then selected for 
comparison. Table 1 shows the number of conformations 
selected for each molecule and the energy range above global 
minimum in which they fall. 

Table 1 



Compound 


Number of CoTiforir.ers 


E.-iers',:-- R.-!-.G" 
Global Minimum 
(Kcal/mol) 


XR905 1 


11 


3.17 


MS-073 


21 


3.49 


GF12091S 


7 


2,68 


CP 1003 56 


n 


5.26 


S97SS 


21 


1.26 


VX-710 


22 


2.99 


B8509-035 


20 


2.93 



"xedminlin" automatically links the programs 
"xedminOl" and "xedf Idprep" . Conformations including those 
5 with XED charges and field points added can be visualized^ in 
"xeddraw" . 

Figure 1 depicts a conformation of XR9051 with XED 

charges added. 

Field points describing the positive, negative and 
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van der Waals surface around the molecule were defined for 
each conformation using the program -xedf Idprep" . Fxgurc 2 
depicts a conformation of XR9051 with related field po.nts. 
IP Figure 2. red represents a positive field point, green 
represents a negative field point and yellow represents a van 
de>- Waals surface field point. 

Comoarxsons between pairs of molecules (21 in total) 
were made, solely on the basis of field point descriptors, 
for each conformation in the sets using the program 
"xedqmfcom". "BCOM" was selected from the menu of 
..xedf^dcreo". This initiate "xedqmfcom". The names of the 
two ^iles (cosmic. fst files) to be compared were entered. 
The results of the comparisons were expressed as coulombic 
overlay energy between field patterns and ranked in energy 
ord-r lo'.'est first using the program "proc.bmf • . The 50 
lowest energv overlays were kept. Table 2 shows the top six 
ranked overlays and related conformational numbers from the 
comparison of XR905i and GF120918. 

Table 2 



Coulombic Overlay 
Enercr.' 



-3945.5 
-3823.7 
-3685.2 
-3639.1 
-3577 6 
-3388 S 



XR9051 
Conformation 
Number 



S 
S 
S 
2 
S 
S 



GF12091S 
Coiitormation 
Number 



5 

6 
3 
1 



Order by 
Coulombic Overlay 



1 

2 

3 

5 
6 



The resc'ts of the comparisons were mapped onto a 
co-^formational mapping diagram and the degree of correlation 
between the field patterns of specific conformations analysed 
uoing the orogram "xedsorf . Figure 3 shows a conformational 
mlpping di'agram describing the top 3 links between molecules. 
Th!s rLals a high correlation between the field patterns or 
specific conformations. No more than^four conformations rrom 
each molecule were sufficient to define the best three 
comparisons between each pair of the seven compounds used. 

The optimum cyclic path between molecules was then 
defined also using the program "xedsort". This facilitated 
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discrimination between conformers of each molecule on the 
optimum cyclic path. Figure 4 depicts the optimum cyclic 
path and related conformations. The conformations that 
appear on the cyclic path were termed "active" conformations 
The "active" conformations were calculated to be 
XR9051 Conformation No. 

GF1.20918 Conformation No. 
MS-073 Conformation No. 

S97G8 Conform.ation No. 

B8509-035 Conformation No 
VX-710 Conformation No. 

CP100356 Conformation No. 
Using the field points of the "active" 
conformations, available from the corresponding cosmic. fst 
files, a pharmacophore was obtained using the "create 
temolate" function within the program "xedsort". The 
parameters used in creating the aggregate average were: 
Electronic hunting distance 3.0 
Surface hunting distance 3.0 
Lov.';3r limit filter 2-0 
The pharmacophore obtained has the following 

narametsrs : 
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.\ton 

Type 

-7 

-7 

-7 

-7 

-7 

-7 

-7 

-7 

-7 

-6 

-6 

-6 

-6 

-6 

-6 



3 
7 



, 381791 
. 379077 
.118141 
. 562676 
. 206078 
-6 .492223 
-6.89920S 
-4 .444944 
-7 .304530 
-6.831137 
5 .201669 
4 .797647 

4 .724869 
-1.074833 

5 .790484 



-0.051772 
-2 .537416 
0.089392 
4 .774324 
1.839680 
-10.028677 
-6.415226 
-7.302401 
8.956919 
-,4.965901 
J). 095870 
5.131311 
0 .054068 
-1.256638 
3.896064 



1.715311 
-5. 803063 
6. 975651 
6. 500367 
-1.274007 
1.675134 
2.683271 
-3.535873 
0.177705 
5.561902 
-2.055756 
-2.954147 
5.399575 
-0.461180 
1.509006 



Charge 

-14 .7214 
-6 . 0644 
-4 . 3289 
-8.6588 
-8.5662 
-2 . 1906 

-7 .7529 

-4 .6301 

-2.2923 

2.3235 

8.2589 

6.4282 
5.1285 

8.2910 

3 .2358 



-6 


-5,261753 


1 .922220 


4 ,363908 


6.6584 


-6 


-1 .473675 


5.651587 


0.238804 


6 . 5755 


-6 


-5,881620 


-0.674191 


-5 . 583631 


3 . 9685 


-6 


0,611350 


-0 .417927 


5 . 020914 


7 . 3589 


-6 


-5 . 783878 


3 .835954 


-4 .612011 


6 . 1148 


-6 


1 .734998 


-0.534123 


10 . 273318 


2 . 7267 


-6 


7 .598676 


-2 .714864 


-1 . 632125 


4 . 5656 


-6 


9 . 587200 


-2 . 109952 


-4 .796533 


At "3 C /t 

4.9354 


-6 


8 .219660 


0.966784 


-8 . 134950 


5.0419 


-6 


4 .675870 


2 .821225 


-7 . 602368 


4 . 9785 


-6 


-1 . 133006 


-9.683890 


0. 147326 


5 . 2926 


-6 


-7. 544412 


-1.565759 


0.572997 


3,8282 


-6 


-9.125527 


-5.337989 


-1 . 040235 


2 . 6320 


-6 


-7.119754 


-8.783302 


-1 .247938 


3.5238 


-6 


2.332763 


-5.381765 


3 .370511 


12 . 2 Db9 


-6 


-6 .407646 


-2 . 912221 


-3 .392133 


3 . 6672 


-6 


-1.476288 


-4 .227778 


1.803737 


7 .4663 


-6 


-9.308277 


4 .502877 


0.694358 


2 . 8187 


-6 


-2 .701459 


-6.564950 


8 . 367678 


8 .4013 


- O 


-0 . 104383 


5.439674 


7.873184 


2 ,2847 


- 5 


0.780038 


0.691121 


1 . 882956 


-11 . 5525 


-5 


-6.232273 


2.601393 


1 .696347 


-10.7761 


-5 


-3 .467500 


-2 .270888 


3 . 913391 


- 1 i . 3 2 3 


-5 


5.578010 


-4 .081225 


0 .238672 


-3 . 5105 


- 5 


-0 . 161253 


0.037560 


-6 .019390 


-8 . 1238 


-5 


■ -2.238659 


-0.779217 


8 .042641 


-8.1531 


-5 


-0 .744474 


-6.033126 


-0.571065 


-8 . 8733 


-5 


-1 .279436 


-5. 982479 


3 . 824223 


-10-0235 


-5 


-7.480469 


-2.520825 


0. 141700 


- 8 . 1691 


-5 


-1,934318 


6.335472 


6 .394043 


-5 • 5024 


-5 


-3 .708762 


5. 991115 


-0.484472 


- 6 , 873 6 


-5 


3.115965 


-2,167775 


-0.311654 


-3 • 9490 


-5 


7 .338803 


2.696229 


-2 . 884882 


- 5 , 4 513 


-5 


- 2,050096 


4.896177 


5. 076291 


-7 . 9992 


-5 


« -4 . 193728 


-2.180720 


-3.354927 


-7 .2592 


-5 


-4 .348258 


1,087337 


-5.092674 


-8.1012 


-5 


-0.521176 


2.639100 


10.414149 


-3-5695 



wherein; atom type describes the nature of the field point, 
-7 represents surface field points, -6 represents positive 



field points and -5 represents negative field points; X, Y 
and Z represent co-ordinate values; and charge represents the 
effective size of the field. Figure 5 is a pictorial 
representation of the pharmacophore. 

In order to show that the pharmacophore produced is 
useful for predicting the MDR modulatory activity of a 
molecule, the field points derived from each of the seven 
known potent MDR modulators used at the outset were compared 
with the pharmacophore. The comparisons were again made 
solely between field point descriptors using the program 

"xedqmfcom" and the results processed using the program 

"proc_bmf " . 

The activities of the MDR modulatory compounds, both 
accumulation IC„ and potentiation IC,,. were then obtained 
either from the literature or by the following assays. 

Assays Fo r Testin g Activity of MDR ?^on n1atory Compounds 

Materials and Methods 
The EMT6 mouse mammary carcinoma cell line and tne 
imR resistant subline AR 1.0 were cultured in RPMI 1640 
0 medium containing 10% foetal calf serum and 2irM glutamine at 
37°C in S-i CO,. Cells were passaged between 1 in 200 and 1 
in 200C in the case of the parental coll lii'-o ana o ^tv.-:. 
in 20 and 1 in 200 in the case of the MDR resistant subline, 
after trypsinisation (0.25% trypsin, 0.2gl\ EDTA) . 
5 n_. Drug accumulation assay 

.\R 1.0 cells were seeded into 96 well opaque culture 
plates (Canberra Packard) . The assay medium contained a 
mixture of tritiated Daunorubicin (DNR). a cytotoxic agent, 
ar.^ unlabelled DNR (0.3 m Ci/ml; 2pM) . The MDR modulators 
0 were serially diluted in assay medium over a range of 

concentrations from 5 nM to 100 ^iM. The cells were incubated 
at 37»C for 1 hr before washing and determination of cell 
associated radioactivity. Results are expressed as % maximum 
accumulation where 100% accumulation is that observed m the 
5 presence of the known RMA verapamil at a concentration of 100 
/iM or as an ICj,,. 

2^ Pot en ti at -inn of Do xorubicin Toxicity 

The MDR modulators were examined for their ability 
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to potentiate the toxicity of doxorubicin in AR 1.0 cells. 
In initial proliferation assays compounds were titrated 
against a fixed concentration of doxorubicin (0.96 //M) which 
alone is non-toxic to AR 1.0 cells. After a four day 
5 incubation with doxorubicin proliferation was measured using 
the colorimetric sulphorhodamine B assay (Skehan et al; J. 
Natl. Cancer Inst. 82 pp 1107-1112 (1990)). 

Compounds which were shown to be able to sensitise 
AR 1.0 cells to 0.86 /xM doxorubicin without high innate 
10 toxicity were selected for further study. Cells were 

cultured for four days with concentrations of doxorubicin 
over the range of 0.01 nM to 50 fxM in the presence of fixed 
concentrations of MDR Modulator at their maximum subtoxic 
levels determined from previous experiments. Proliferation 
15 was quantified as described by Slcehan et ai, loc cit. The 

IC50 (concentration required to reduce proliferation to 50% of 
the untreated controls) for doxorubicin alone and for the MDR 
Modulators were derived and used to calculate the 
potentiation index (PI) : 
20 IC50 for Doxorubicin alone 

PI= 

IC50 for Doxorubicin plus modulator 



Table 3 





Potentiation 
IC50 nM 


Data 
Source 


Compound 


Accumulation 
IC50 nM 


GF120918 


5-33 


8-20 


1.2 


S9788 


32 


1580 


1 


CP100356 


160 


300 


2.3 


MS-073 


1000 


-200 


4.3 


B8509-035 


200 


1000 


1 


VX-710 


2500 


-2000 


2 


XR9051 


180 


200 


1 
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Table 4 
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Example No 


EHD 


SHD 


LLF 


(r) 

Log [acc IC50] 
overlay energy 


(r) 

Log [pot IC50] 
vs Coulombic ! 
overlay energy 


2 


2.0 


0.5 


2.0 


0 .54 


0.63 1 


1 3 

r 


2.0 


3.0 


2.0 


0 .79 


i 

0.85 1 


i 


3.0 


0.5 


2.0 


0 .54 


0.63 j 




3.0 


1.0 


2.0 


0 .54 


0.63 j 




2.0 


0.5 


5.0 


0 .56 


0.65 1 




2.0 

i 


1.0 


5 . 0 


0 .55 


i 

0.69 ( 


i ' 


2 . 0 


2.0 


5.0 


0 . 72 


0.71 I 




3.0 


0.5 


5.0 


0 . 56 


0.65 


(1 10 
1 


3 . 0 


1.0 


5.0 


0.59 


0.69 i 


t , . 


3.0 


3.0 


5.0 


0.77 


0.80 



1.. 
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r T. A I M S 



1. A pharmacophore comprising a 3 -dimensional array of 
field points defining a shape and volume wherein the set of 
3-dimensional field points is the aggregate average of the 

5 field points derived from a plurality of molecules which 
possess multi drug resistance (MDR) modulatory activity. 

2. A pharmacophore according to claim 1, wherein the 
field point? are derived by: 

(i) subjecting each molecule to conformational hunting 
10 to generate a plurality of conformers for each molecule; 

(ii) applying extended electron density (XED) charges to 
each molecular conformer generated in (i) ; 

(iii) defining field points describing the positive, 
negative and van der Waals surface around each conforr;ier; 

15 (iv) making an intermolecular energy comparison between 

the field points defined in (iii) ; 

(v) obtaining the three lowest energy comparisons for 
each intermolecular comparison and mapping them on a 
conformational mapping diagram,- 

20 (vi) defining an optimum cyclic path which relates the 

molecules; and 

(vii) defining the field points for the conformer of each 
molecule v/hj.ch lies on the optimum cyclic path defined in 

(vi) . 

25 3. A pharmacophore according to claim 1 or 2 wherein 

the plurality of molecules which possess r4DR modulatory 
activity includes one or more of XR9051, MS-073, GF120918, 
CP100356, S9788, VX-710 and B8509-035. 

4 . A pharmacophore according to any one of claims 1 to 

30 3, which when compared by an overlaying m.ethod with the field 
poiiits of each of the plurality of compounds from which the 
pharmacophore is derived gives coulombic overlay energy 
values which when plotted against log [accumulation IC50] has 
a sample correlation coefficient (r) of not less than 0.5. 
35 5. A pharmacophore according to any one of claims 1 to 

4 having the parameters 
Atom 

Type X Y Z Charge 

-7 -1.381791 -0.051772 1.715311 -14.7214 
40 -7 4.879077 -2.537416 -5.803063 -6.0644 

-7 -4.418141 0.089392 6.975651 -4.3289 



I- V ' 



1- - 



-7 3.562676 

-7 7,206078 

-7 -6.492228 • 

-7 -6.899206 

-7 -4,444944 

-7 -^7,304530 

-6 -6.831137 

-6 5.201669 

-6 4.797647 

-6 4.724869 

-6 -1.074833 

-6 5.790484 

-6 -5.261753 

-6 -1.473675 

-6 -5.881620 

-6 0.611350 

.6 •5.7S3878 

-6 1.734998 

-6 7.59867G 

-6 9.537200 

-G 8.219660 

-6 4.675870 

-6 -1.133006 

-6 -7.54»W12 

-6 -9.125527 

-6 -7.119754 

-6 2.332763 

-6 -6 . 407646 

- 6 - 1 . 476283 

-S -9.308277 

-6 -2.701459 

-6 -C.104383 

-5 0.780038 

-5 -6.232273 

-5 -3.467500 

-5 5,578010 

-5 -0.161253 

-5 -2.238659 

-5 -0.744474 

-5 -1.279436 

-5 -7.480469 
-5 -1.934318 
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4,774324 
1.839680 - 
•10.028677 
•6.415226 
-7.302401 
8.956919 
-4 .965901 
0.095870 
5.131311 
0 . 054068 
-1.256638 
3 . 896064 
1 . 922220 
5.651587 
-0 . 674191 
-0 . 417927 
3 .835954 
-0.534123 
-2.714864 
-2 . 109952 
0.966784 
2 .821225 
-9.683890 
-1 .565759 
-5 .337989 
-8 .783302 
-5 . 381765 
-2.912221 
-4 .227778 

4 . 502877 
-6 . 564950 

5 .439674 
0 .691121 
2 .601393 

-2 ,270888 
-4 .081225 

0,037560 
-0-779217 
-6.033126 
-5.982479 
-2.520825 

6.335472 



6.500367 
-1,274007 

1 ,675134 

2 .683271 
-3.535873 

0 . 177705 

5,561902 
-2.055756 
-2 . 954147 

5.399575 
-0 .461180 

1 . 509006 

4 .363908 
0 .238804 

-5 . 583631 

5 . 020914 
-4 .612012 

10 .273318 
-1,632125 
-4 .796533 
-8 . 13495C 
-7 . 802366 

0 . 147326 
0 . 572997 

-1 . 040235 
-1-. 247938 
3 . 370511 
-3.392133 

1 . 803737 
0 .694358 
8 , 367678 
7.873184 
1 .882956 
1.696347 
3 . 913391 
0 .238672 

-6 .019390 
8 .042641 

-0.571065 
3.824223 
0 . 141700 
6.394043 



-8.6588 
■8.5662 
-2,1906 
-7.7529 
-4.6301 
-2 .2923 
2.3235 
8 .2589 
6.4282 
5.1285 
8 .2910 
3.2358 
6 . 6584 
6.5755 
3.9685 
7.3589 
6.1148 
2 .7267 
4 .5656 
4.9354 
5.0419 
4 . 9785 . 
5.2926 
3.8282 
2.63z0 
3.5238 
12.2559 
3 .6672 

7 .4663 
2 .8187 

8 .4013 
2 .2847 

-11,5525 

-10 .7761 

-11.2323 

-3 . 5105 

-8.1238 

-8 .1531 

-8.8733 

-10.0235 

-8.1691 

-5.5024 



28 

-5 -3.708762 5.991115 -0.484472 -6.8736 
-5 3.115965 -2.167775 -0.311654 -3.9490 
-5 7.338803 2.696229 -2.884882 -5.4513 
-5 2.050096 4.896177 5.076291 -7.9992 
5 -5 -4.193728 -2.180720 -3.354927 -7.2592 

-5 -4.348258 1.087337 -5.092674 -8.1012 
-5 -0.521176 2.639100 10.414149 -3.5695 

wherein; atom type describes the nature of the field point. 

-7 represents surface field points, -6 represents positive 
10 field points and -5 represents negative field points; X, Y 

and Z represent co-ordinate values; and charge represents the 

effective size of the field. 

6. A method of screening a molecule for MDR modulatory 

activity, which method comprises; 
15 (i) obtaining field points describing the positive, 

negative and van der Waals surface around the molecule; 
(ii) measuring the energy of overlay when the molecular 

field points obtained in (i) are fitted into a pharmacophore 
as defined in any one of claims 1 to 5; and 
20 (iii) selecting the molecule as a candidate MDR modulator 
if the energy of overlay measured in (ii) is equal to or 
lower than the energy of overlay between field points for 
molecules known to possess MDR modulatory activity and the 
said pharmacophore. 
25 7. A method of designing a molecule having MDR 

modulatory activity, which method comprises the step of 
determining the structure of a compound which has field 
points that have, when fitted into a pharmacophore as defined 
in any one of claims 1 to 5. an energy of overlay equal to or 
30 lower than a predetermined value. 

B A method according to claim 7 wherein the 

predetermined value is the energy of overlay measured when 
the field points for a plurality of molecules kno-...i to 
possess MDR modulatory activity are fitted into the said 

35 pharmacophore . . 

9 A compound which possesses MDR modulatory activity ^ 

identified using a pharmacophore as defined in any one of - 
claims 1 to 5 or by a method as defined in claim 6 or 
designed using a method as defined in claim 7 or 8. 
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